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Abstract 

In two dimensions a Hopf bifurcation occurs as a Spiral Point switches from stable 
to unstable (or vice versa) and a periodic solution appears. There are, however, more 
details to the story than this: The fact that a critical point switches from stable 
to unstable spiral (or vice versa) alone does not guarantee that a periodic 
solution will arise, ^ though one almost always does. Here we will explore these 
questions in some detail, using the method of multiple scales to find precise conditions 
for a limit cycle to occur and to calculate its size. We will use a second order scalar 
equation to illustrate the situation, but the results and methods are quite general and 
easy to generalize to any number of dimensions and general dynamical systems. 



^Extra conditions have to be satisfied. For example, in the damped pendulum equation: x+^x + sinx = 0, 
there are no periodic solutions for ! 
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1 Hopf bifurcation for second order scalar equations. 



1.1 Reduction of general phase plane case to second order scalar. 

We will consider here equations of the form 

X + h{x,x, fi) = , (1.1) 
where h is a smooth and is a parameter. 

Note 1 There is not much loss of generality in studying an equation like (1-1), as 
opposed to a phase plane general system. For let: 

x = f{x,y,iJ.) and y = g{x,y, fi) . (1.2) 

Then we have 

x = fxx + fyy = fxf + fyg = F{x, y, ii) . (1.3) 
Now, from x = f{x, y, fi) we can, at least in principle,'^ write 

y = G{x,x,iJ.) . (1.4) 

Substituting then (1-4) into (1-3) we get an equation of the form (1-1).^ 

1.2 Equilibrium solution and linearization. 

Consider now an equilibrium solution"^ for (1.1), that is: 

x = X{ii) such that h{0, X, //) = , (1.5) 

^We can do this in a neighborhood of any point {x^,y^) (say.a critical point) such that fy{x^,y^, fi) ^ 0, 

as follows from the Implicit Function theorem. If = 0, but gx 7^ 0, then the same ideas yield an equation 

of the form y + h{y,y,fi) = for some h. The approach will fail only if both fy = = 0. But, for a 

critical point this last situation implies that the eigenvalues are fx and gy, that is: both real ! Since we are 

interested in studying the behavior of phase plane systems near a non— degenerate critical point switching 

from stable to unstable spiral behavior, this cannot happen. 

^Vice versa, if we have an equation of the form (1.1), then defining y hy y = G{x,x,fj,), for any G such 

that the equation can be solved to yield x = f{x, y, fi) (for example: G = x), then y = G^x + GxX = g(x, y) 

upon replacing x = f and x = —h. 

H.e.: a critical point. 



18.385 MIT (Rosales) Hopf Bifurcations. 



4 



so that X = X is a solution for any fixed fi. There is no loss of generality 



in assuming 
(1.6) 



X{fi) =0 for all values of fi 

since we can always change variables as follows: Xq[^ = X{p.) + Xnew- 

The linearized equation near the equilibrium solution x = (that is, the equation for x 
infinitesimal) is now: 

X — 2ax + f3x = , (1.7) 

where 



a = a{ii) = -i/ii(0,0,/i) 



and 



(3 = (3{fi) = h,{0,0,fi) 



The critical point is a spiral point if /5 > a^. The eigenvalues and linearized solution 
are then 



(where co = \/ ^ — a^) and 



where a and to are constants. 



A = q; ± zcj 



X = ae" cos {uj{t — to)) 



(1.8) 
(1.9) 



1.3 Assumptions on the linear eigenvalues needed for a Hopf bi- 
furcation. 

Assume now : At /i = the critical point changes from a stable to an unstable spiral 

point (if the change occurs for some other fi = fi^, one can always redefine /^qU = /ic + A'new)- 
Thus 

[a < for < and a > for > 0, with /? > for small. 
In fact, assume ; 

• I. h is smooth. 



II. ck(0) = 0, ;9(0)>0 and ^a(0)>0.s 



(1.10) 



We point out that, in addition, there are some restrictions on the behavior of 
the nonlinear terms near the critical point that are needed for a Hopf bifurcation 
to occur. See equation (1.22). 



^This last is known as the Transversality condition. It guarantees that the eigenvalues cross the imaginary 
axis as /i varies. 
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1.4 Weakly Nonlinear things and expansion of the equation near 
equilibrium. 

Our objective is to study what happens near the critical point, for small. Since for = the 
critical point is a linear center, the nonlinear terms will be important in this study. Since 
we will be considering the region near the critical point, the nonlinearity will be weak. 



Thus 



we will use the methods introduced in the Weakly Nonlinear Things notes. | 



For x^x, and // small we can expand h in (1.1). This yields 



X 



+ ulx + {\Ax^ + Bxx + \Cx'^] + 

+ \ {Dx^ + '^Ex'^x + '^Fxx'^ + Gx^} 

- 2p^xii + Qx/i + 0{e'^, e^/i, eji^) = , 



(1.11) 



where we have used that h{0, 0, /i) = and a{0) = 0. In this equation we have: 
d 

A. = _/j(o,0,0) = /3{0) > 0, with > 0, 



B. A: 

C. 

D. 



dx^ 

1 



dxdx 

d 



E. e is a measure of the size of {x,x). Further: both e and fi are small. 



1.5 Explanation of the idea behind the calculation. 

We now want to study the solutions of (l-H). The idea is, again: for e and // small the 
solutions are going to be dominated by the center in the linearized equation x + UqX = 0, 
with a slow drift in the amplitude and small changes to the period^ caused by the higher 
order terms. Thus we will use an approximation for the solution like the ones in section 2.1 
of the Weakly Nonlinear Things notes. 



^We will not model these period changes here. See section 2.3 of the Weakly Nonlinear Things notes for 
how to do so. 
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1.6 Calculation of the limit cycle size. 

An important point to be answered is ; I What is epsilon? | (1.12) 

This is a parameter that does not appear in (1.1) or, equivalently, (1-11). In fact, the only 
parameter in the equation is // (assumed small as we are close to the bifurcation point fi = 0). 
Thus: 



I f must 1)6 related to (1-13) 



In fact, f will be a measure of the size of the limit cycle, which is a property of the 

equation (and thus a. function of // and not ai'l)itrai'y all). 
However; 



We do not know e a priori! How do we go aljout determining it? 



I The idea is; | If we choose e "too small" in our scaling of {x,x), then we will be looking 
"too close" to the critical point and thus will find only spiral-like behavior, with no limit 
cycle at all. Thus, we must choose e just large enough so that the terms involving 
in (1.11) (specifically 2p^^x, which is the leading order term in producing the sta- 
ble/unstable spiral behavior) are "balanced" by the nonlinearity in such a fashion that 
a limit cycle is allowed. In the context of Two-Timing this means we want fi to "kick 
in" the damping/amplification term 2p^^x at "just the right level" in the sequence of 
solvability conditions the method produces. Thus, going back to (1-11), we see that^ 

• The linear leading order terms x + u)qX appear at 0(e). 

• The first nonlinear terms (quadratic) appear at O(e^). 

However; Quadratic terms produce no resonances, since sin^ ^ = -(1 — cos 29) and 
there are no sine or cosine terms. The same applies to cos^ Q and to 
sin Q cos Q. 

• Thus, the first resonances will occur when the cubic terms in x play a role =^ we must 
have the balance 

0{x^) = 0{fix) , (1.14) 



This is a crucial argument that must be well understood. Else things look like a bunch of miracles! 
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1.7 The Two Timing expansion up to O(e^). 

We are now ready to start. The expansion to use in (1.11) is 

X = exi{T,T) + e'^X2{T,T) + €^X3{t,T) + . . . , (1.15) 
where < e ^ 1, 27r-periodicity in T is required, T = coot, coq is as in (1.11)^, r is a slow 



time variable and e is related to fi by 



where | = ±1 ] (which u we take depends 



on which "side" of = we want to investigate). 

What exactly is r? Well, we need r to resolve resonances, which will not occur until the cubic 



terms kick in into the expansion =^ \ t = eH'. (This is exactly the same argument used to 
get (1.14)). 

(1.11) becomes: 



Then, with 



d 

dT 



lOqx" + lOqX + ||A(Jo(a:')^ + Bujqxx' + |Ca:^| + 

\{Dujl{x'f^?,Eujl{x'fx^2,Fujax'x''^Gx^ + (1.16) 
2e2(jo4 - 2t^vp^u)Qx' + e^vilx + O(e^) = . 
The rest is now a computational nightmare, but it is fairly straightforward. Without 
getting into any of the messy algebra, this is what will happen: 



At 0(e) + = 0. Thus 



Xi = ai(T)e* + c.c. 



(1.17) 



for some complex valued function aiij). We use complex notation, as in the Weakly Non- 
linear Things notes. 



At 0{^) 



'^o {^2 + 2:2} + {quadratic terms in Xi and x\} = Q . 



(1.18) 



From the first bracket in (1.16), the quadratic terms here have the form: 

where Ci and C2 are constants that can be computed in terms of Wq, ^, B and C. 

Since the solution and equation are real valued, C2 is real. Here, as usual, * indicates the 
complex conjugate. 



^Same as the linear (at /i = 0) frequency. No attempt is made in this expansion to include higher order 
nonlinear corrections to the frequency. 
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No resonances occur and we have 



(1.19) 



At 0(e3) 



(Jo (xg + xg) + 2a;oa:'i^ - 2up^coox\ + z^J^Xi + CNLT = . 



(1.20) 



where CNLT stands for Cubic Non Linear Terms, involving products of the form X2X1, 
x^Xi, X2x[, x'2x[, (x'l)^, {x[)'^Xi, x\x\ and x\. These will produce a term of the form 
da\a{e^'^ -\- c.c. plus other terms whose T dependencies are: 1, e^^*^ and e^^*^, none of 



which is resonant (forces a non periodic response in ^3). Here 



c? is a constant that can be computed in terms of Wq, A, B, C, D, E, F and G . I (1-21) 



This is a big and messy calculation, but it involves only sweat. In general, of course, 
lm{d) 7^ 0. The case lm{d) = is very particular, as it requires h in equation (l.l)to be just 
right, so that the particular combination of its derivatives at a: = 0, i; = and = that 
yields Im(d) just happens to vanish. Thus 



Assume a nondegenerate case : lm{d) ^ 



(1.22) 



For equation (1.20) to have solutions x^ periodic in T, the forcing terms proportional to e^*^ 
must vanish. This leads to the equation: 



2iOoi—ai — 2vp^ijjQiai + vQ.ai -\- d a\ ai = 0. 
dr 



(1.23) 



Then write 



This yields 



and 



ai = pe*^ , with p and 9 real , p > . 



■^p = vp^{l - vqp^)p., 



(1.24) 



(1.25) 



where 
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Equation(1.24) provides a correction to the phase of xi, since xi = 2pcos(r + 9). The 
first term on the right of (1.24) corresponds to the changes in the linear part of the phase 
due to // 7^ 0, away from the phase T = loqI at // = 0. The second term accounts for the 
nonlinear effects. 

The second equation (1.25) above is more interesting. First of all, it reconfirms that for 
< (that is, = — 1) the critical point (p = 0) is a stable spiral, and that for > (that 
is, = 1) it is an unstable spiral. Further 



If Im((i) > 



If lm{d) < 0. 



Then a stable limit cycle exists for 

> (i.e. u = 1) with p = ^j2uJop'^{lm{d))-^ . 
Supercritical (Soft) Hopf Bifurcation. 

Then an unstable limit cycle exists for 
li<0 (i.e. u = -1) with p = ^ -2iOop'^ {lm{d))-^ 
Subcritical (Hard) Hopf Bifurcation. 



(1.26) 



Notice that p here is equal to — the radius of the limit cycle. 



1.7.1 Remark on the situation at the critical bifurcation value. 

Notice that, for p = (critical value of the bifurcation parameter)^ we can do a two timing 
analysis as above to verify what the nonlinear terms do to the center. The calculations are 
exactly as the ones leading to equations (1.23)-(1.25), except that = and e is now a small 
parameter (unrelated to /i, as = now) simply measuring the strength of the nonlinearity 
near the critical point. Then we get for p = — radius of orbit around the critical point 



2e 



(1.27) 



From this the behavior near the critical point follows. 

^Then the critical point is a center in the linearized regime. 

^°This is the way one would normally go about deciding if a linear center is actually a spiral point and 
what stability it has. 
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\m{d) > <^=^ Soft bifurcation <^=^ Nonlinear terms stabilize. 

For [1 = critical point is a stable spiral. 

\m{d) < <;=^ Hard bifurcation <^=^ Nonlinear terms de-stabilize. 

For 11 = critical point is an unstable stable spiral. 



Clearly < 



1.7.2 Remark on higher orders and two timing validity limits. 

As pointed out in the Weakly Nonlinear Things notes, Two Timing is generally valid for some 
"limited" range in time, here probably |t| <^ e~^. This is because we have no mechanism 
for incorporating the higher order corrections to the period the nonlinearity produces. If 
we are only interested in calculating the limit cycle in a Hopf bifurcation (not it's stability 
characteristics), we can always do so using the Poincare-Lindsteadt Method. In particular, 
then we can get the period to as high an order as wanted. 

1.7.3 Remark on the problem when the nonlinearity is degenerate. 

What about the degenerate case 



lm{d) = 



In this case there may be a limit cycle, or there may not be one. To decide the question 
one must look at the effects of nonlinearities higher than cubic (going beyond 0{€'^) in the 
expansion) and see if they stabilize or destabilize. If a limit cycle exists, then its size will 
not be given by •y'^l/ij, but something else entirely different (given by the appropriate balance 
between nonlinearity and the linear damping/amplification produced by a 7^ when ^ 
in equation (1.7)). The details of the calculation needed in a case like this can be quite hairy. 
One must use methods like the ones in Section 2.3 of the Weakly Nonlinear Things notes 
because: even though the nonlinearity may require a high order before it decides the issue of 
stability, modifications to the frequency of oscillation will occur at lower orders. We will 
not get into this sort of stuff here. 



^^Note that Re{d) ^ in (1.24) produces such a change, even if Im(d) = and there are no nonlinear 
effects in (1.25). 



